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Turning attention to the fatigue process itself, Elber (ref. 1 1) observed that fatigue-crack surfaces 
contact each other even dining tension-tension cyclic loading. This simple observation of the crack-closure 
phenomenon immediately explained many crack growth characteristics that had been troubling researchers 
during the ’60s. Since then, several closure mechanisms have been identified, among them, "plasticity" 
induced closure. These new closure mechanisms and the influence of the plastic wake on the local crack-tip 
strain field have greatly advanced an understanding of the fatigue-crack growth and fracture behaviour of 
metallic materials, under, for example, variable amplitude loading. Furthermore, the occurrence of crack- 
closure significantly affects the local crack driving force stress-intensity factor and plays a crucial role in 
quantifying the fatigue crack growth or arrest characteristics of the material in question. Thus, while the 
qualitative interpretation of the closure phenomena justifies a large number of fatigue crack anomalies, a 
quantitative knowledge of the crack closure stress intensity level is required to correlate fatigue crack 
growth data. 

Most fatigue design requirements based on damage tolerance concepts assume the existence of flaws 
in the component either from the initial delivery or from some later stage during the service life. To assure 
a high level of reliability and before these flaws grow to critical lengths, it is necessary to either repair or 
replace the component. Thus, in-service inspections are required to detect various sizes and shapes of 
cracks and other defects. For a reliability analysis of an in-service airframe, for example, fatigue, 
environment and accidental events are three sources of damage that must be taken into consideration. 

Hence, the importance of fatigue; the unifying nature attributed to closure; the need for a stochastic 
process interpretation of the whole fatigue process, and the need for reliability analyses are clear. While 
crack closure and stochastic processes are briefly detailed, this paper primarily concerns an investigation of 
new methods of assessing the fatigue reliability of structures based upon probabilistic approaches. It 
summarizes; i) the experimental determination of a statistically significant number of crack growth rates for 

a Ti-6A1-4V titanium alloy; ii) the determination of the Closure-Lognormal, c, m, <7 2 , £ parameters for Ti- 
6A1-4V; iii) the simulation of the crack propagation based upon these Closure-Lognormal parameters; iv) 
the determination of the McGill-Markov A and k parameters for the Ti-6A1-4V alloy in question; and v) the 
RAM assessment of Ti-4A1-6V. 


Fatigue Crack Closure 


An understanding of crack closure under cyclic tension was developed by Elber (ref. 11) who 
showed that the occurrence of premature contact between the opposing crack faces during unloading was 
due to the residual plastic stretch in the crack wake. During the loading portion of a cycle, the elastic 
constraints acting on the residual material in the wake of an advancing crack, keep the crack tip closed until 
these constraints are overcome by the externally applied load. The stress intensity factor associated with a 
fully opened crack, K op , based upon the crack opening load, P op , is necessary for the quantitative 

knowledge of the effective stress intensity range factor, A K eff , which is used in life prediction through 
K eff = K nmx - K op and the well known Paris-Erdogan relation da / dN = cAK eff m . 

A K eff is an appropriate field parameter for correlating crack growth under constant-amplitude 

loading conditions to the influence of a large number of variables, such as an actual load spectrum, load 
ratio and/or a specific environment, that are known to affect the rate at which cracks grow in any practical 
situation. As an example, fatigue crack closure effects at different load ratios have been extensively 
investigated by Ritchie et al. (ref. 12), who initially showed a significant difference in the da / dN - A K 
behaviour due to different load ratios over a wide range of A Ks, but, when plotted as a function of A K eff , 
were able to eliminate this discrepancy and consolidate the curves into a narrow band. 
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Stochastic Processes 


multitude of posable outcomes in such a way that statistical regularity can be oteerS 
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The McGill-Markov Process. In order to further evaluate the stochastic characteristics of fatigue 
crack propagation processes, it is necessary to consider the concepts behind the McGill-Markov process. 
Any Markov process is based upon the assumption that the prediction of the process that is going to occur 
is influenced only by the properties of the current state which the process is in and not by the history that 
led to its present state. According to Provan and Rodriguez (ref. 16) and the nature of fatigue crack 
growth, the crack propagation process is a discrete-state continuous-parameter and nonhomogeneous 
Markov process since the transition probability density is a variable which depends only on the time 
difference r. Furthermore, the crack size variable, a(t), can only be measured to within equipment and 
operator limitations. By considering the observable zones, i, the crack size may be written as: 

a, < a{t) < a, + Aa ; , ( 5 ) 

where i is the state number and Ao i is the width of a state. A discrete- state and continuous-parameter 
stochastic process, e t}, may be described by a one step transition probability of the form: 

P{a(t) = j\a(T) = i} = p ij (T,t), 0 <T<t, (6) 

where i and j are integer states and T and t are times. This probability is called the transition probability 
and is defined as the probability of a transition from state i to state j during the time interval % to t. 

In order to solve the Kolmogorov differential equations which govern p i]t an infinitesimal transition 

scheme must be specified. After a review of the existing intensity functions. Provan and Rodriquez (ref. 
16) developed a new intensity which gives a good description of the time evolution of material property 
degradation processes. The intensity functions of the McGill-Markov process are: 
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X and k are empirical system parameters which are determined by a fit to experimental data. They are 
positive empirical system parameters that describe all of the various effects that influence the process, such 
as, temperature, material properties, experimental error, etc. When these parameters are found, the system 
may be modeled. If the system is changed, new system parameters must be found. These intensity 
functions are then used in the Kolmogorov forward differential equation, which becomes: 


dpjjif, 0 

dt 
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Hence, by solving this Kolmogorov equation governing a linear nonhomogeneous birth process, the 
transition probabilities p, y (T,i) are determined. 


Finally, in order to determine the entire history of the crack propagation distribution, the total 
probability may be continuously monitored via the fundamental absolute probability relation: 


P jit) = 

;=i 


( 10 ) 
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where P (t) is the probability of being in the j-th (later) state at time t; />(r) is the probability of being in 

for anv ^ ^ 113118111011 Probabilities of a Markov process. Hence, 
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furthermore there will be only one inspection/correction process in a given interval time. Hence, the 
question: "what is the optimum time for this procedure ?" is an appropriate one. In order to decrease the 
probability of failure an inspection/correction procedure will be earned out at some time. An inspection too 
early in the service life will, on the one hand, remove few components that may subsequently fail while a 
later inspection may be too late to remove components that will have failed. The optimum time for 
inspection will depend on several variables such as: critical crack size, repair size, inspection process, and 
the quality of replacement components. Hence, the McGill-Markov model, used in conjunction with a 
failure control methodology, can be a useful tool for obtaining valuable reliability information. 


EXPERIMENTAL PROCEDURE AND BASIC RESULTS 


Exerimental Procedure. A computer controlled, increasing stress intensity factor A K test method, 
according to the American Standard Test Method ASTM E647 (ref. 18), was applied to eighteen standard 
C(T) specimens manufactured from a forged Ti-6A1-4V jet engine, fan disk grade, titanium alloy. The 
main objectives of carrying out these experiments were to examine the stochastic properties of crack growth 
for Ti-6A1-4V, and then to analyze the results in an effort to establish the parameters associated with the 
Closure-Lognormal and McGill-Markov stochastic processes and in the reliability assessment of 
components manufactured from this titanium alloy. 


The fatigue tests were carried out under the control of the in-house "FATIG" computer program 
based on a crack closure compliance method. For the compliance calculations, the CMOD measurement 
was determined using a mechanical clip gage with a maximum resolution of 0.00025 mm, while that of the 
load cell was 0.005 kN. The FATIG program utilizes three main loops and determines: i) the crack length 
based on load vs. CMOD data (compliance); ii) the A K based upon the ASTM E647 standard; and iii) the 
A K efI based upon A K and the closure load. During each "data acquisition" block, the load vs. CMOD 
curve for 200 individual data points was obtained. The lower limit was a variable such that the nonlinearity 
was distinguished by FATIG as the closure load. A linear curve fit was made through the remaining points 
to obtain the normalized compliance, and the relationship between the compliance and crack length was 
determined by that proposed by Mirzaei and Provan (ref. 19), namely: 


EBv 

P 


3(ct + l)(a + 2.3) 
a 2 - 2a + 1 


+ 8(1 + v)a; 



( 12 ) 


In eqn(12), a is the normalized crack length, E = 117,000 MPa is the elastic modulus of Ti-Al-4V, v is 
the CMOD, B the specimen thickness, P the load and W the specimen width. From this information the 
A K eff was determined from the expression given in the ASTM E674 standard, except that A P = P mn — P op , 

P np being the closure load. 


When a specific test was completed, the stored data was analyzed with the second in-house "FADA" 
program to obtain a da / dN vs. A K eff curve using the incremental polynomial method with n=3 , i.e., 7 
successive data points. The FADA program plotted the crack length, a, vs. the number of cycles, N, and 
the crack growth rate, da / dN vs. AK eff . 


Experimental Test Results For Ti-6Al-4V. Results were obtained from eighteen standard C(T), Ti- 
6A1-4V specimens prepared according to conventional procedures. The fatigue crack propagation for all 
eighteen specimens are shown in fig. 1(a), while the associated da / dN vs. A K eff curves are presented in 
fig. 1(b). Examination of the data revealed a few surprising observations. The most unexpected 
observation was the change in the growth rate which occurred in almost every test, indicating that fatigue 
crack propagation is not a stable, smooth, well ordered process. A close observation of the fatigue crack 
growth rates show two transition points, namely, at A K eff = 8 — 9 MPa and at A K eff s 13 — 14 MP a . 
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ence, the fatigue crack growth process is highly complex, especially in the Paris-Ergodan regime and it is 

from thp le - t0 ha y e t d f e rentPans-Erdoga n regimes correlated to these experimental results. However, 

the first ™ir nt ^ abllit y anal Y ses the ^ rst of the fatigue crack growth rate curve, i.e., prior to 
the first transition point, is more important than the second. * 


THE CLOSURE-LOGNORMAL CHARACTERISTICS OF Ti-6A1-4V. 


As introduced above, a lognormal process is defined by the fact that its logarithm is a normal 

Stefs^ function tfthe Sf ge ? C J al w ° f fatigue . crack Propagation laws indicates that the crack growth 
and o n l,! 4 slt y, fac tor, maximum stress intensity factor, stress amplitude, load ratio, 

and so on Some commonly used crack growth rate functions, such as the Paris-Erdogan model (ref. 20 s ) 
re such that c and m are functions of the load condition and environment. 

? ss £ ssm ent Of the closure phenomena necessitates the evaluation of the state of 
residual stress and strain m the neighbourhood of the crack tip and the extent of crack closure. Crack 

closure effects are most pronounced at low A K levels. This can be seen in fig. 2 which shows the results 
o tamed from the eighteen specimens of Ti-6A1-4V tested as described above. Hence, P o for Ti-6A1-4V 
as a function of crack length, a, is found by a polynomial curve fit to be as follows : 


P op = 2. 90203- 0140572a + 0.0021723a 2 ; 


(kN;mm). 


(13) 


lr . t A S ^fl ai !^? re y. iousl J; the crack closure concept may be used to describe the influence of the actual 
load spectrum, load ratio and/or environmental parameters and, therefore, the quantitative knowledge of the 

"closu!^^ ? f Vel t 1S r T kCd t0 COrrelate fatigue CTack growth rate da ta. Using the § sense of 

closure , m the form of the effective stress intensity range factor for fatigue crack propagation the 

^ qn i 14 ^- 1S ado P ted m such a manner that it is independent of component geometry, 
oadmg spectrum or load ratio, environment, etc. In this way, the entire description of the statistical scatter 
m fatigue crack growth data may be incorporated into the assessment of the reliability of any prospective 
component prior to its manufacture. y y P ros P ectlve 

suggested* 106 ' m ° rd6r t0 aCCOunt for * e random nature o f the crack growth rate, the following model is 


~ = X(t)c(AK a y 


(14) 


where a(t) is now the random crack size and X(t) is a non-negative random process. In this case, A K 
logStoofXVS!rf^n/l°4 U k Mows Aa d t X(t) * " ‘° gn0m,al random process ' By ^ the 


is a 


log (— ) = logX(0 + logc + mlog(AK eff ). 


(15) 


By substitution T = log(^), Z(r) = logX(r), C = Iogc, r=log(A/r„), eq„(15) may be 
rearranged into the form: 


= mT + C + Z(t). 


(16) 
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Z(t) describes the inherent scatter of a specific material. These parameters, as well as the variance d 2 and 
the autocorrelation parameter, £ , are obtained from the Ti-6A1-4V test results of the crack growth rate vs. 
the effective stress intensity range in the logarithmic scale. 


According to the stochastic model, eqn(14), and the sense of closure for fatigue crack propagation, 

four important Closure-Lognormal parameters, namely m, c cx 2 and £ , may be interpreted as material 
properties. Table 1 shows these Closure-Lognormal parameters for Ti-6A1-4V. Hence, a complete 
specification of a material’s stochastic fatigue crack growth characteristics are found from the Closure- 
Lognormal model. These parameters may be interpreted as material properties that are independent of 
component geometry, loading spectrum or load ratio, environment, etc. In this way, the entire description 
of the statistical scatter in actual test results or simulated fatigue crack growth data (see the next section), 
may be incorporated into the assessment of component reliability. Since, by using these parameters and the 
simulation procedure, a definitive reliability, availability and maintainability procedure may then be earned 
out using the associated McGill-Markov parameters, a complete specification of a material s stochastic 
crack growth characteristics based upon the Closure-Lognormal interpretation of scatter is of paramount 
importance to the assessment of crack growth rates in any practical situation. 


SIMULATION OF CRACK PROPAGATION - Ti-6A1-4V 

In order to remove the limitation of a specific crack geometry, loading, or environment, etc., and to 
be in a position to use both the previous Closure-Lognormal and the following McGill-Markov procedures 
of describing the statistics associated with the fatigue crack growth process, a simulation procedure plays a 
crucial role in transferring information from data generated by a standardized Closure-Lognormal procedure 
to real crack propagation processes as they occur in actual components and under any of a large number of 
specified situations. The stationary Gaussian random process Z(t), may be simulated using the well-known 

Fast Fourier Transform (FFT) technique to simulate the two parameters, £ , and a 2 . 

As an example, the C(T) specimen geometry along with a knowledge of the P op for the Ti-6A1-4V 

was used to generate the simulated fatigue crack growth characteristics of this material. The results of these 
analyses are shown in fig. 3, which clearly shows that the two sets of data, one experimentally obtained 
and the other simulated, superimpose on each other. 

A K eff , incorporating a quantitative knowledge of the crack opening stress level, now appears as an 
appropriate field parameter for correlating knowledge of the constant-amplitude crack rates to practical 
growth rate situations. Specifically, knowing in terms of m, c d 2 and £ , the crack growth rate as a 
function of A K eff for a specific material, such as Ti-6A1-4V, then the influence of an actual load spectrum, 

load ratio and/or environment may be incorporated during the design stage into an estimate of the fatigue 
reliability of a specific component. In this way, the simulation process plays an important role in removing 
the dependency on actual experimental results obtained under specific loading and environmental 
conditions. Hence, the closure behaviour, either experimentally determined or predicted by a model, 
expressed in terms of the effective stress intensity factor range expressions found in handbooks and the 
Closure-Lognormal parameters, constitute the only information required for analysing any particular 
situation. 


THE MCGILL-MARKOV PARAMETERS FOR Ti-6 A1-4V 

The major objective of this section is the determination of the McGill-Markov stochastic fatigue 
crack growth properties of Ti-6A1-4V. As was described in the introduction, the two constant parameters 
A and k in the intensity functions < 7 ; and q ijy defined in eqns(7) and (8), are related to the stochastic 
properties of the material system being investigated. Accordingly, the probability of the crack tip being in 
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howev^S^S h^t terS if? dct ermi'Kd by a fit to the experimental or simulated data. Several steDS 
data m ^ be ? re * ,s lteratlve P roc ^ » undertaken. The first step is the normahzatt P on of 

data to an mitral crack length of a, at hme t=0. This is done to eliminate the crack initiation stage. In the 

case, a Siam srae of 0.4mm along with a failure state of N, = 40 were chosen tewasjudg^ tot a 

-mevaluesof A and .which give 


X = 0.14 


JC = 0.98. 


(17) 


fig 4 showrEth? mT* e f Pf rimen tal,simulated and McGill-Markov model predictions, shown in 

rlSVnAtTr T* who , le concept of treating fatigue crack growth as stochastic processes is a flexihS 
X and ^ PI C ng CTaCk P ™ pagation for “dustiially significant materials and situations. The fact that the 

iv^innmE am< dT f a ? p P bc * b Jf ln reliability prediction and inspection/correction procedure 
development, as detailed m the following section, again allude to the benefits of this a^pS 


RELIABILITY ANALYSIS OF Ti-6A1-4V 
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nr „ cpn th< L results for a change in acceptable reliability level from 0 9999 to 0 9995 are 
(state) from N r = 20, through N r = 25, to = 30 prior to 2.0E+06 cycles. There are a total of two 
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inspection/collections for N, = 20, three for N, = 25 and four for N, = 30. These figures illustrate the 
type of information that can be obtained from this approach. 


Inspection Optimization. Finally, in order to carry out an inspection optimization analysis, the 
design engineer must control such variables as failure state, N r repair size (state), N r , and the desired 

level of reliability. By using N f = 40 and N r = 30, the solid line in fig. 8 is obtained. From this figure it 

is apparent that the optimum time for inspection/correction is at 1 .700E+06 cycles and that the total 
probability of failure is decrease by 88% over the no inspection case. By changing repair policy, such that 
N is varied from 30 to 20, different curves, also shown in fig. 8, are obtained. In this way, it is shown 
that while the optimum inspection time is increased by increasing the repair crack size from 20 to 30, the 
overall reliability of system is reduced. 


CONCLUDING REMARKS 


The conversion of AKinto A K eff through the inclusion of closure effects plays a sigmficant role m 

reliability analysis. By using AK e# , it now becomes possible to transfer the stochastic properties of crack 

growth rates, measured under ideal laboratory conditions, to practical situations '^ 

closure into both the fatigue crack growth rate description and the lognormal interpretation of the scatter has 

led to di^devetopment of the Closure-Lognormal model which describes the statistical nature of crack 

growth rates. 

On the other hand, the McGill-Markov process, employing data generated by simulations of the 
information contained in the Closure-Lognormal interpretation of the basic material s 
characteristics, vary with respect to closure effects that describe the influence of variations in the loadmg, 
environment, crack geometry, etc. With the crack propagation characteristics being predictable, the 
reliability and inspection processes may then be evaluated. 
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Figure 1(a): a vs N for Ti-6A1-4V 


Figure 1(b): da/dN vs A K tff forTi-6Al-4V 








Crack length (mm) 


Figure 2: Closure Load vs Crack Length for Ti-6A1-4V 


Table 1: Closure-Lognormal Parameters for Ti-6A1-4V 
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Figure 5: Reliability as a Function of Fatigue Cycles 



a) At 0.9999 maintenance reliability 



b) At 0.9995 maintenance reliability 
Figure 6: Inspection/Correction Procedure for Ti-6A1-4V 
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Figure 7: Inspection/Correction under different Repair Sizes 
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Figure 8: Inspection Optimization for Ti-Al-4V 



